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Various systems of very different physical nature exhibit interesting singular be- 
havior in the vicinity of their critical points. Their correlation functions reveal self- 
similar behavior with universal critical dimensions: they depend only on few global 
characteristics of the system (like symmetry or space dimensionality). Quantitative 
description of critical behavior is provided by the renormalization group (RG). In the 
RG approach, possible types of critical regimes (universality classes) are associated 
with infrared (IR) attractive fixed points of renormalizable field theoretic models. 
Most typical phase transitions belong to the universality class of the O n -symmetric 
if) 4 model of an n-component scalar order parameter. Universal characteristics of the 
critical behavior depend only on n and the space dimensionality d and can be calcu- 
lated in the form of the expansion in e = 4 — d or within other systematic perturbation 
schemes; see the monograph [T] and the literature cited therein. 

It has long been realized that the behavior of a real system near its critical point 
is extremely sensitive to external disturbances, geometry of the experimental setup, 
gravity, presence of impurities and so on; see e.g. the monograph [2] for the general 
discussion and the references. What is more, some disturbances (randomly distributed 
impurities or turbulent mixing) can produce completely new types of critical behavior 
with rich and rather exotic properties, like e.g. expansion in y/e rather than in e; see 
e.g. 0. 

These issues become even more important for the nonequilibrium phase transitions, 
because the ideal conditions of a "pure" stationary critical state can hardly be achieved 
in real chemical or biological systems, and the effects of various disturbances can 
never be completely excluded. In particular, intrinsic turbulence effects can hardly 
be avoided in chemical catalytic reactions or forest fires. One can also speculate that 
atmospheric turbulence can play important role for the spreading of an infectious 
disease by flying insects or birds. Effects of different kinds of turbulent and laminar 
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motion on the critical behavior were studied e.g. in [3] [9]. 

In this paper we study effects of strongly anisotropic turbulent mixing on the crit- 
ical behavior of two paradigmatic models: the equilibrium model A, which describes 
purely relaxational dynamics of a nonconserved scalar order parameter (see e.g. [1]), 
and the Gribov model, which describes the nonequilibrium phase transition between 
the absorbing and fluctuating states in a react ion- diffusion system [10]. The velocity 
is modelled by the <i-dimensional generalization of the random shear flow introduced 
in [XT] within the context of passive scalar advection. 

In the Langevin formulation the models are defined by stochastic differential equa- 
tions for the order parameter if) = if>(t, x): 

d t ip = A {(-r + d 2 )ifj - V{iP)} + c = 0, (1) 

where dt = d/dt, d 2 is the Laplace operator, A > is the kinematic (diffusion) 
coefficient and r oc (T — T c ) is the deviation of the temperature (or its analog) from 
the critical value. The nonlinearity has the form V(if>) = uif> 3 /3\ for the model A and 
Vfy) = 5 l '0 2 /2 for the Gribov process; g and u > being the coupling constants. The 
Gaussian random noise ( = x) with zero mean is specified by the pair correlation 
function: 

(C(t,x)C(t',x')> = 2A5(t-t')5 (d) (x-x') ( 2 ) 

for the model A and 

(C(t,x)C(t',x')> =^A^,x)5(t-^ (d) (x-x') (3) 

for the Gribov process. The factor if) in front of the correlator (CD) guarantees that in 
the absorbing state the fluctuations cease entirely, while the factor 2A in ([2]) ensures 
the correspondence to the static ip 4 model. 

Coupling with the velocity field v = {vi(t, x)} is introduced by the replacement 
d t — > V( = d t + Vid{, where Vt is the Lagrangian (Galilean covariant) derivative. 

Let n be a unit constant vector that determines distinguished direction ( "direction 
of the flow"). Then any vector can be decomposed into the components perpendicular 
and parallel to the flow, for example, x = xj_ + nx» with xi ■ n = 0. The velocity 
field will be taken in the form v = m;(t,x_|_), where w(t,x_|_) is a scalar function 
independent of x». Then the incompressibility condition <9jt>j = d\\v(t, xj_) = is 
automatically satisfied. 

For v(t, x_|_) we assume a Gaussian distribution with zero mean and the pair cor- 
relation function of the form: 

(v(t,x ± )v(t',x' ± )) = D5(t-t>) J ^^exp{ik ± -(x ± -x'J}-^. (4) 
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Here k± = [kj_ [ , D > is a constant amplitude factor and £ an arbitrary exponent, 
which (along with the conventional e = 4 — d) will play the part of a formal RG 
expansion parameter. The IR regularization is provided by the cutoff k± > m in <^j. 
The natural interval for the exponent is < £ < 2 with the most realistic Kolmogorov 
value £ = 4/3; the "Batchelor limit" £ — > 2 corresponds to smooth velocity. 

According to the general theorem (see e.g. [TJ), stochastic problems described 
above can be reformulated as field theoretic models of extended set of fields $ = 
The role of the coupling constants is played by the parameters 

w,# 2 ~A 4 - d , w = D/\~At, (5) 

where A is some typical ultraviolet (UV) momentum scale. From (JSJ) it follows that 
the both models become logarithmic (all the coupling constants are simultaneously 
dimensionless) at d — 4 and £ = 0. Thus the UV divergences in the Green functions 
manifest themselves as poles in e = 4 — d, £ and, in general, their linear combina- 
tions. The divergences can be removed by the standard renormalization procedure. 
In order to ensure multiplicative renormalizability of the model, it is necessary to 
split the Laplacian in (TTJ into the parallel and perpendicular parts d 2 — > dj_ + fd 2 
by introducing a new parameter / > (in the anisotropic case, these two terms are 
renormalized in a different way). Then the corresponding RG equations are derived; 
their IR attractive fixed points determine the IR asymptotic scaling regimes of the 
models. Detailed analysis of the model A can be found in [8] and for the Gribov 
model it will be given elsewhere. Below we only present the results of the explicit 
one- loop calculation (leading order in £ and e), which appear similar for the both 
models. There are four fixed points: 

1) Gaussian (free) fixed point: g* = = u>* = 0, IR attractive for e < 0, £ < 0. 

2) g* = = (exact result to all orders), u>* ~ £; IR attractive for £ > and 
£ > 2s. In this regime, the nonlinearity V(ip) in the stochastic equation ([T]) becomes 
irrelevant, and we arrive at the model of a linear convection-diffusion equation for a 
passive scalar field ip. 

3) = (exact), g\ ~ ~ e; IR attractive for e > 0, £ < for the model A 
and £ > 0, £ < e/12 for the Gribov model. In this regime, effects of the turbulent 
mixing are irrelevant, the isotropy violated by the velocity ensemble is restored and 
the leading terms of the IR behavior coincide exactly with those of the standard 
models 

4) The most interesting fixed point with nontrivial positive values of g*, u* and 
w*. It is IR attractive for e > 0, £ > for the model A and e > 0, £ > e/12 for 
the Gribov model. It corresponds to new nontrivial IR scaling regimes, in which both 
nonlinearities in the stochastic equations and the turbulent mixing are important; 
the corresponding critical dimensions depend essentially on the both RG expansion 
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parameters e and £ and are calculated as double series in these parameters. This 
behavior reveals strong anisotropy and belongs to new, completely nonequilibrium, 
universality classes. The realistic values £ = 4/3 and d = 2 or 3 belong to these 
classes. 

The regions of IR stability for these fixed points are shown in Fig. 1. For com- 
parison, the boundaries for the isotropic stirring [9] are given. In the one-loop ap- 
proximation, all the boundaries are given by straight lines; there are neither gaps nor 
overlaps between the different regions. The latter fact is valid to all orders of the RG 
expansion, although the boundaries between the regions 3 and 4 (determined by the 
nonlinearity V(ip)) and 2 and 4 (determined by the velocity statistics) will become 
curved. 

Let us illustrate the consequences of the RG analysis for the spreading of a cloud 
of the "agent" in the turbulent environment. The mean-square displacement Ri(t) in 
the i-th direction at time t > of the agent's particle, which started from the origin 
x' = at time if — 0, is given by the relation 

ti*(t)= f dxa$G(t,x), G(t,x) = (V>(£, x)^ f (0, 0)), x=|x|. (6) 

The linear response function G(t, x) has the following asymptotic representation 

G(t,x) = x'^'^F (x ± r x / A ", x\$-^I*»,tx£) , (7) 

where x± = |xj_| and F is the scaling function. The set of critical dimensions A* for the 
fields and parameters is determined by the fixed point. In particular, for the Gribov 
model at the fixed point 4 in the one- loop approximation one has A w = 2 — (2e — £)/23, 
A,, = 1 + (12£ - e)/23, A^ = A^ t = d/2 + (14f - 5e)/46, A T = 2 + 3(f - 2e)/23. 
Substituting © into (jSJ) at r = (that is, directly at criticality) gives the power laws 

R\{t) <xt a ±, Rl(t)<xn, (8) 

where the exponents a±, an are simply related to the dimensions A*. For the fixed 
point 3 (linear passive advection) one obtains exact results a± = 1, a\\ = 1 + £/2. For 
the transverse direction this gives the ordinary diffusion law R±(t) oc t 1//2 , while in the 
direction of the flow the spreading is accelerated and for £ = 4/3 takes on the form 
R\\(t) oc t 5//3 , or equivalently dR^/dt oc R^/ 5 ■ This "4/5" law differs from the classical 
"4/3" Richardson's law for the isotropic case. For the regime 4 the exponents in (jSJ) 
depend on the both parameters £ and e. 
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£ = 2e (Shear flow) 



£ = e/2 (Isotropic flow) 



£ = e/12 (Gribov model) 
£ = (Model A) 



Figure 1: Regions of stability of the fixed points in various models. 
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